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y ! Lifetimes of the first two metastable states in Sc^+ and are determined usins; 

the relativistic coupled-cluster theory. There is a considerable interest in studying 

, the electron correlation effects in these ions as though their electronic configurations 

, are similar to the neutral alkali atoms, their structures are very different from the 

^SJ ' latter. We have made a comparative study of the correlation trends between the 

CN ■ above doubly ionized systems with their corresponding neutral and singly ionized 

. iso-electronic systems. The lifetimes of the excited states of these ions are very 

o ; 

important in the field of astrophysics, especially for the study of post-main sequence 

o ■ 

' evolution of the cool giant stars. 

$H ' PACS numbers: 



I. INTRODUCTION 

With the recent progress in high performance computational techniques, it is now possible 
to carry out accurate calculations on heavy atomic systems using powerful relativistic many- 
body theories . It is not possible to treat all the heavy atomic systems by a single 
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many-body theory due to the differences in their electronic structures. The relativistic 
coupled-cluster (RCC) theory has been successfully applied earlier to a wide range of atomic 
properties J] including the lifetime calculations of Pb"*" jsl and alkaline-earth ions 6]. The 
role of electron correlation effects is dramatic in some of those cases. In the present work, 
we have carried out the lifetime calculations of the lowest (^5/2 and Si/2 states in Sc III and 
Y III ions. 
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FIG. 1: Schematic diagram of the energy levels and transition lines between the ms and nd states 
in Sc^^ and Y^^. The n and m represent the principal quantum numbers of the ground and 
metastable states where n = 3 and m = 4 for Sc^+ and n = 4 and m = 5 for Y^+ . 



Generally, the forbidden transition lines find numerous interesting applications in astro- 
physics which includes the determination of the elemental abundances in different celestial 
objects and also in the laboratory astrophysics and plasma physics studies. The spectral 
lines of many alkaline-earth ions including Y and its doubly charged ions have been observed 
in the spectra of many hot stars like Ap stars [7] and these ions seem to be over abundant 
in Hg-Mn stars compared to their normal stellar abundances. This anomaly observed over 
abundances than that expected from the theoretical estimations based on the s-process of 
neutron capture is studied extensively in the literature j^, ol. Thus, the study of adjacent 
elements of the same row of the periodic table like, the Sr-Y-Zr triad of elements, can be 
used to understand the nucleo-synthetic origin of these elements and it can reveal the post- 
main-sequence evolution of the chemically peculiar stars, such as Hg-Mn stars fz^. The Y 
.11 and Sc III ions can also be used for the diagnostics and modeling of the stellar plasma 
lo| . Recently, Brage et. al. 11| have used the multi-configuration Hartree-Fock (MCHF) 
method in calculating the lifetimes of the 5p, 5d and 6s states of Y III which were measured 



earlier by Maniak et. al. [12]. However the lifetimes of the 4^5/2 and 5si/2 states in Y 
III seem to be unreported so far. Though there are a large number of results available for 
the transition wavelengths of Sc 111, it seems that there are no measurements on either the 
lifetimes or the transition probabilities of the 8^5/2 and 4si/2 states. For this reason, we have 
carried out the calculations of the lifetimes of the lowest ^5/2 and Si/2 metastable states in 
Sc 111 and Y 111 employing a powerful many-body theory known as the RCC theory. As the 
experimental wavelengths are known to high accuracy we have used them in our calculations 
in order to minimize the errors in our lifetime estimations. The low-lying excited states and 
their electromagnetic decay channels are shown in the figure [H Even though Sc 111 and 
Y III belong to K- and Rb- iso-electronic sequence, respectively, the ground states of the 
former two are ^^3/2 states (n = 3 for Sc III and n = 4 for Y III) unlike the latter two 
for which the ground states are {n + 1) si/2- Again though the considered systems are the 
doubly charged ions, they have the same ground states as that of their neutral counterparts. 
Thus the ions considered in this work are interesting from the theoretical and as well as 
astrophysical point of view. We analyze the role of correlation effects in these calculations 
and compare their trends with the previously studied lifetimes of the (^5/2 states [l^ in alkali 
atoms and alkaline earth metal ions [6]. 

II. THEORY AND METHOD OF CALCULATION 

The transition probabilities for the spontaneous emissions due to E2 and Ml electromag- 
netic transitions from the state i to the state / are given by, 

,P2 1.11995 X lO^s^po 

Ml _ 2.69735 X 10^3 

where, [Jj] = 2Jj + 1 is the degeneracy of the state i, S = \{f\0\i)\'^ is the transition line 
strength in atomic unit (au) for the corresponding transition operator O and wavelength A 
(A). 

As shown in the figure [T|, the 6/5/2 state can decay to the ^3/2 state in the considered 
systems via E2 and Ml channels; hence the total transition probability for the state 
can be expressed as, 

And5/2 = ^^d5/2^nd3/2 + ^nd5/2^nd3/2- (2-3) 
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Similarly, the S1/2 state can decay to the 6^3/2 state in these systems via E2 and Ml channels 
and decay to the ^5/2 state via E2 channel; the resulting total transition probability is 

A — _|_ 4E2 I aE2 fc) A\ 

^msl/2 — ^msl/2~^nd3/2 ' ^msl/2^nd3/2 ' ^msl/2^nd5/2- K^^-^J 

In the above expressions, n and m represent the principal quantum numbers of the d and 
s states; where, n = 3 and m = 4 for Sc 111 and n = 4 and m = 5 for Y 111. The single 
particle reduced matrix elements of E2 and Ml operators used in our calculations can be 
found from p^] . The lifetimes of these states are determined by taking the inverse of their 
total probabilities. 

We use the following RCC ansatz to construct the wave functions for one valence states 
in the considered systems; 

\^,) = e^{l + S,}\<l>,), (2.5) 

with the reference state 1$^,) = aJ,|$o) for the valence electron v and I'I'o) is the Dirac-Fock 
(DF) state of the closed-shell system constructed by the Dirac-Coulomb (DC) Hamiltonian. 
This is the relativistic method of the universal valence coupled-cluster theory proposed by 
M...e,ee e. a. ^ a.d ,a.e. co.pacU. fo^ulated by L,nd..en Q. ,„ .Ke above RCC 
wave function expression, we define T and as the closed-shell and open-shell excitation 
operators which excite core electrons and valence with core electrons, respectively. We have 
considered all single, double and a subset of important triple excitations approximation in 
our calculations which is known as the CCSD(T) method in the literature 0, It has 
been found in many studies in the similar systems as considered here that the CCSD(T) 
method is powerful enough to account the correlation effects to evaluate various properties 



to high accuracy j^, ^ ll7|, ll8|]. This is primarily because of the fact that this approach 
considers the most important correlation effects like core-polarization and pair-correlation 
effects to all-orders jl8| . 

The transition matrix elements for the operator O are calculated in this approach using 
the following expression; 



fi 



{<^f\{l + S\]e''Oe^{l + S,]\^,) 
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{<i>f\{i + s}}o{i + sm) 

, (2.6) 



where, O = e^''Oe^. Since the above expression contains many non-hnear terms those 
will be computationally expensive if they are calculated directly, we follow the extended 
Wick's theorem [l9| to break these terms in simple form as discussed below. Although these 
procedures are given in detail earlier [s], [3], we mention them briefly here. First, we divide 
O as the effective one-body, two-body etc. operators and later these terms are sandwiched 
between the necessary operators. To save the computational time, we calculate and save 
only the one-body terms which are later used where other terms are computed directly in 
the final property calculations. It has to be noticed that the fully contracted O terms do not 
contribute in the present calculations. Contributions from the normalization of the wave 
functions (Norm.) are accounted by. 

Norm. = (^.101 $i){^^=-l}, (2.7) 

where, 

= (<l>o|e^'e^|<l>o) + (<f „|{5]e^'e^^4|<l>„) 
for the corresponding valence electrons v (= i, /). 

III. RESULTS AND DISCUSSIONS 

In Table [U we present the wavelengths and transition line strengths of the low-lying 
transitions from the lowest ^5/2 and Si/2 states to the ground states of Sc III and Y III. 
The transition wavelengths for Sc III are taken from the critically evaluated atomic data 
base of National Institute for the Standard and Technology (NIST) j2o| and the transition 
wavelengths for Y III are taken from 21| for the given transitions. The calculated transition 
strengths are used to determine the corresponding transition probabilities which are given 
in the last column of the above mentioned table. 

The calculated total transition probabilities for the lowest ^5/2 and Si/2 states are given 
in Table [TTl The lifetimes are determined by taking the reciprocal of these quantities for the 
corresponding states. The Ml transition probabilities between the fine structure levels, i.e. 
the dr,/2 — > (^3/2 transitions are 5-6 orders of magnitude larger than the E2 transitions, where 



as, the Ml transition probabilities for the 6^3/2 —>■ Si/2 transitions are 5-6 orders smaller than 



6 



TABLE I: Transition wavelengths, line strengths and transition probabilities (As) of different low- 
lying transitions in Sc^'^ and Y^"*". 



Transition 

(f ^i) 




A 

(in A) 


Transition 
channel 


Strength 
f au) 


A 

sec~^ 


Sc2+ 


















Ml 


9 3751 


R '?A'iA F,-5 








E2 


2.7164 


1.5291 E-11 




3d 2Z)3/2 


3916 


Ml 


7.9495E-7 


1.786 E-4 








E2 


12.9157 


7.8584 


As 25i/2 - 




3946 


E2 


19.5007 


11.413 


Y2+ 












Ad 2l)5/2 - 


> 4d 2Z)3/2 


138093 


Ml 


2.3822 


4.0667 E-3 








E2 


9.9573 


3.7011 E-8 


5s 2^1/2 ^ 


4d 27^3/2 


13392 


Ml 


1.1136E-7 


6.2544 E-7 








E2 


37.6554 


4.895 E-2 


5s 25i/2 ^ 


4d 2^5/2 


14830 


E2 


57.7644 


4.5090 E-2 



the corresponding E2 transitions (cf. Table [T]) for both the ions. Thus, the Ml transitions 
play crucial roles in determining the lifetimes of the ^5/2 states but it is the other way around 
for the Si/2 states as the E2 transitions play crucial roles in their lifetime determinations. 
As it can be seen from Table HTl that the lifetime of the 8^5/2 state in Sc III is very large 
implying that it is a highly forbidden state. However, the lifetime of the Asi/2 state is a 
few fraction of a second. On the other hand, the lifetime of the 4^5/2 state in Y III is 246 
seconds which again is a long lived metastable state, where as, the 5si/2 state has a lifetime 
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TABLE II: The transition probabilities and the hfetimes of the Si/2 and ^5/2 states in Sc^ and 
Y2+. 



State 


Transition Probabihty 


Lifetime 




(sec""'") 


(sec) 


Sc2+ 






3d 2^5/2 


8.2434 E-5 


12130.86 


45 25i/2 


19.2715 


0.0519 


Y2+ 






Ad 21^5/2 


4.0668 E-3 


245.89 


5S 2Si/2 


9.4041 E-2 


10.63 



of only 10 seconds. This also indicates that the correlation effects in the calculations of Ml 
transition amplitudes of the ^5/2 —>■ transitions are very small and hence these results 
can be ascribed as relativistic in origin. The transition probabilities of the states given in 
Table HTl will be very useful in astrophysics for the accurate determination of the abundances 
of these elements. We estimate that our lifetime results can have an error bar of less than 
2%. The source of errors includes the truncation of the coupled-cluster series at the double 
excitations and partial inclusion of triple excitations, the omission of Breit and QED effects 
in our calculations. 

To understand the importance of the correlation effects, we present the contributions 
from the individual terms of Eq. (12.61) in Table IIIII to Table |Vl The contribution given by 
the operator O is the Dirac-Fock term and the difference between O and O represents the 
pure core- valence (CV) correlation effects. The remaining large contributing terms OSu and 
082% with their conjugate terms represent the all-order pair-correlation and core-polarization 
effects, respectively, for the given valence electron. 

The reduced matrix elements of the Ml and E2 transition operators for the — > (^3/2 
transitions in Sc2+ and Y2+ are given in Table IIIII The transition amplitudes due to E2 are 
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TABLE III: Reduced matrix elements of the Ml and E2 transition operators from ^5/2 states to 
the ground states in au in the considered systems. 



RCC 




Sc2+ 






Y2+ 




terms 


Ml 




E2 


Ml 




E2 



Contributions from effective one-body terms of O. 






-1.55 


-1.93 


-1.55 


-3.54 


0-0 


1.78E-2 


1.37E-2 


1.30E-2 


1.98E-2 




5.41E-4 


9.50E-2 


7.05E-4 


1.16E-1 


Sip 


3.16E-4 


9.47E-2 


-8.33E-5 


1.13E-1 


OS2i 


-3.44E-5 


4.34E-2 


-3.35E-4 


6.99E-2 


sip 


1.88E-5 


4.23E-2 


2.06E-4 


6.98E-2 


s\pSu 


-3.5E-3 


-1.14E-2 


-2.11E-3 


-1.15E-2 


S\pS2i 


1.41E-6 


3.34E-4 


-1.28E-6 


5.57E-4 


S^^OSxi 


-9.73E-7 


3.02E-4 


1.74E-7 


5.08E-4 


SlpS2i 


-3.33E-2 


-1.92E-2 


-2.85E-2 


-3.86E-2 



Contributions from effective two-body terms of O. 



TlOS2i 1.57E-3 9.66E-4 1.38E-3 1.73E-3 

SlpTi 1.57E-3 9.65E-4 1.37E-3 1.72E-3 

Tl0S2i 0.00 0.00 5.68E-7 -1.25E-5 

slpT2 0.00 0.00 -6.29E-7 -1.21E-5 

Norm. 2.28E-2 2.44E-2 1.99E-2 4.06E-2 

Total -1.54 -1.65 -1.54 -3.16 
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TABLE IV: Reduced matrix elements of the Ml and E2 transition operators from Si/2 states to 
the ground states in au in the considered systems. 



RCC 




Sc2+ 






Y2+ 




terms 


Ml 




E2 


Ml 




E2 



Contributions from effective one-body terms of O. 



o 


3.75E-7 


4.05 


-2.56E-6 


-6.70 


o-o 


-8.77E-4 


-6.70E-3 


l.OlE-4 


4.85E-3 


OSu 


-2.61E-5 


-3.91E-1 


7.64E-6 


4.00E-1 


sip 


-2.17E-5 


-2.68E-2 


-7.74E-6 


7.83E-2 


OS2i 


2.26E-6 


-3.25E-2 


-1.06E-5 


5.87E-2 


sip 


1.72E-6 


5.74E-3 


-1.71E-5 


1.92E-2 


s\pSu 


-7.34E-7 


l.llE-2 


-9.91E-8 


-1.93E-2 


S\pS2i 


3.59E-8 


-8.21E-4 


-2.92E-7 


1.49E-3 


sl^OS\i 


2.93E-8 


1.90E-3 


-1.07E-7 


-2.00E-3 


SlpS2i 


1.90E-5 


2.64E-2 


-4.09E-4 


-5.75E-2 



Contributions from effective two-body terms of O. 



TlOS2i 1.50E-8 -7.09E-4 -3.58E-8 1.60E-3 

SlpTi -8.90E-9 -7.74E-4 1.55E-8 1.63E-3 

Tl0S2i 0.00 0.00 -1.38E-9 4.89E-6 

slpT2 -3.01E-7 3.45E-3 2.24E-7 -4.79E-3 

Norm. l.llE-5 -4.46E-2 4.22E-6 7.76E-2 

Total -8.92E-4 3.59 -3.34E-4 -6.14 
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TABLE V: Reduced matrix elements of the E2 transition operator for the corresponding si/2 to 
(^5/2 states in au in the considered systems. 



RCC 


Sc2+ 


Y2+ 


terms 


E2 


E2 



Contributions from effective one-body terms of O. 



O 

o-o 
slfO 

OS2i 

sip 

S\pSu 

S\pS2i 

SyOSli 

slfOS2i 



4.97 
-8.91E-3 
-4.80E-1 
-3.31E-2 
-4.01E-2 
8.80E-3 
1.36E-2 
-l.OlE-3 
2.36E-3 
3.37E-2 



-8.28 
7.10E-3 
4.86E-1 
9.86E-2 
7.25E-2 
1.75E-2 
-2.36E-2 
1.82E-3 
-2.50E-3 
-7.25E-2 



Contributions from effective two-body terms of O. 



Tl0S2i -8.71E-4 1.96E-3 

slpTi -9.50E-4 2.00E-3 

Tl0S2i 0.00 5.59E-6 

SlpT2 4.23E-3 -5.86E-3 

Norm. -5.46E-2 9.51E-2 



Total 



4.42 



-7.60 



11 



slightly larger than the amplitudes due to Ml transitions, but the latter totally decide the 
lifetimes of the ^5/2 states (see Table [T]). It is because, the transition probabilities due to 
E2 transitions depend on fifth power of their wavelengths in the denominator, where as, the 
Ml transition probabilities depend on third power of their wavelengths in the denominator. 
The total contribution of the correlation effects, which can be evaluated as the difference 
between the DF and the CCSD(T) results, in the calculation of Ml transition amplitudes 
are around 0.5% for both the systems, however, they are more than 10% in the case of the 
E2 transition amplitudes. 

The reduced matrix elements of the Ml and E2 transition operators of the transition, 
•Si/2 in Sc^"^ and Y^"^ are given in Table HVl The transition amplitudes due to E2 are 

almost 3-4 orders of magnitude larger than the Ml amplitudes and the former give largest 
contributions to the lifetimes of the Si/2 states. The correlation effects are very large in both 
the Ml and E2 transitions in Sc^^, where as, they are large only in the case of Ml transition 
in Y2+. 

In Table IVl we present the E2 transition amplitudes of the si/2 — > ^5/2 transitions. The 
trend of the contributions from different RCC terms are similar to the ^5/2 — > Si/2 transitions 
of the corresponding iso-electronic alkali metal atoms [13|] and alkaline earth metal ions j^, 
but these amplitudes are smaller in the considered systems. The correlation effects in these 
calculations are just around 12% and 9% in Sc III and Y III, respectively. Again, the 
correlation effects in Ml transition amplitudes of the —>■ transitions are very small 
compared to the 6^3/2 6pi/2 transition amplitude in Pb"*" [18\, where we observe similar 
behavior for the E2 transition amplitudes. 

The determination of the lifetimes of the Si/2 states depend crucially on the accurate 
calculation of the E2 transition amplitudes of both si/2 (^3/2 and Si/2 —>■ ^5/2 transitions, 
where as, the Ml amplitudes of the former transitions play negligible roles. 



IV. CONCLUSION 

We have determined the lifetimes of the lowest ^5/2 and Si/2 metastable states in Sc^+ 
and Y^+ which seem to be unavailable in the literature. The transition strengths and 
the transition probabilities of these systems find interesting applications in astrophysics 
in understanding the post-main sequence evolution of the cool giant stars. These doubly 
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ionized ions may also be suitable for other ion trapping experiments as they are reported to 
have long lifetimes, which need to be explored by the experimentalists. 
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